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rr^YtTg^T-Trnreatrr^ - - — 

Tho  rooulto  «'£  a  two-yoor  program  on  dynamic  fracture  bohnvlor  of  rocka  are 
doocrlbed,  The  goal  was  to  develop  tho  capability  to  predict  t:Iie  fragment  olzo  dlstri| 
button  of  rock  resulting  from  known  dynamic  IobcIs,  This  goal  hoa  boon  largely 
realized  for  n  simple,  well-characterized  rock  type,  Arkansas  novaculitc,  under  one- 
dltaon/tionnl-strain  impact  loading,  A  gaa  gun  wss  used  to  accolorcte  flat  plates 
against  flat  rock  specitieno,  Ytterbium  plezoresistl vo  stress  gages  were  used  to 
tjoaauro  stress  hlstorlco,  end  tho  dynatslc  tensile  r.trongth  and  fragment  slzo  dtetrl- 
butions  were  determined.  Tho  mochanlea  of  fragmentation  was  deduced  from  f rcetogruphl£ 
obocrvntloas  on  impacted  cpocloens. 

A  computational  f ragnontotlon  model  was  developed  that  treats  quantitatively  tho 
lour  stages  of  tho  hypothesized  fragmentation  mechanism:  (1)  activation  of  inherent 
flaws,  (2)  crack  growth,  (3)  crack  coalcsconco,  and  (4)  fragmentation.  Thin  jsodcil  was 
luovi'lcd  into  a  ono-dimenslonni ,  finite  difference  wavo  propagation  computer  code  to 
obtain  o  capability  to  predict  tho  fragment  size  distribution. 

The  required  input  parameters  Include  the  load  history  r.nd  such  reck  specific 
properties  os  tho  Initial  flaw  size  distribution,  tho  fracture  toughness,  and  the 
crock  growth  velocity.  Those  material  parameters  were  determined  for  Arkansas 
noveculite,  and  n  calculation  vas  tando  to  vinulnto  tho  conditions  ol  r,  dynamic  impact 
experiment.  The  calculated  and  experimental  fragment  size  distributions  (Figure  12) 
are  In  qualitative  agreement,  end  Indicate  that  it  13  feasible  to  make  successful 
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Tho  results  ©2  o  swo-yoat?  prorata  ©a  dynstle  fraetUFe  behavior  of 

rmko  aro  defleribcd.  'tho  gesi  woo  to  dovalop  the  capability  t<*  predict 

the  fFcgisent  ela©  distribution  of rsefc  resulting  trm  Stnasrs  dynselc  loads. 

•  *  .  '  ! 

This  goal  has  boon  largely  reBliced  for  n'clsrplo,  tseil-charisct orised  task 

typo,  Arfecnseo  Rovaenlitc,  under  ©ne-diK«Bs4©aol«ctraia  lespect  lottOies* 

h  0^*  B«sa  was  osKifi  to't<eeci©»ete:  f?-at  plates '  o^r-ilof:?.  f'taz  ©salt  npseVstns. 

Vtterblusa  pioz^oelstive  ©trees  gag©©  wos5©  used  i©  coeoare  o  trout?  historic^ 

and  the  dynvnie  tonal lo  screegth  end  irogoeat  oico  distributions  were 

dot ©rained. ;  ■  Tho  ocstenlcta  ©2.  fragcentcfciea-  was  deduced  frea  irsote^ropbio 

©boevvnuemc  e»  iapoctcd  opecineas. 

A  computational  f ragoantulioa  mi'"5cl  tvn»  developed  thrst  trooto  quanti¬ 
tatively  tho  four  stages  of  tho  hypothooised  Iroocontation  eactiunion: 

(1)  netivation  of  inherent  flaw,  (2)  eraeh  growth,  <3)  crock  cooXonecneo, 
end  (4)  frograentotion,  Thin  eodel  won  inserted  inS©  o  one-diEonaSesn?,, 
finite  difference  wove  propoyetlon  compute*.  code  to  ootain  a  capability 
to  predict  the  fraguoiit  eiso  distribution. 

Tho  required  input  paronoters  ineludo  tho  load  history  end  ouch  *>©ch 
specific  propertioo  ao  the  initial  flow  size  distribution,  tho  fracture 
toughness,  and  the  crack  growth  velocity.  These  material  pnroasotoro  were 
determined  for  Arkansas  novaculito,  ond  o  calculation  van  made  to  ciraulnte 
tho  conditions  of  a  dynamic  impact  exporiraont.  Tho  calculated  end  experi¬ 
mental  fragment  size  distributions  (Figure  12)  arc  In  qualitative  agree¬ 
ment,  and  indicate  that  It  in  feasible  to  make  successful  quantitative 
predictions  of  rock  fragmentation  based  on  knowlodge  of  a  few  measurable 
rock  proportion.  Charectorl zntion  work  on  Sioux  quartzite,  llolcton 
limestone,  and  Wontorly  granite  ‘voe  begun. 
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1  )  NTftOOliCT  I  ON 

Rock  fragmentation  under  high  rate  load  applications  la  o  lit.tl*- 
understood  phenomenon .  Thera,  is  *t  present  no  satisfactory  theoretical 
heals  lor  predicting  dynamic  .fracture  behavior,  although  auch  a  capability 
would  be  vary  crnful  in  the  eolutloh  of  many  practical  mining:  aiul  civil 
engineering  problot!. s .  For  example ,  rapid  excavation  could  be  done  more 
safely  and  ocoflosticslly,  tho  stability  of  structure*  In  rook  could  bo 
designed  and  evaluated  with  morn  confidence,  and  the  efficiency  of  rock 
disintegration  processes  could  be  improved.  It  was  with  thl*  Motivation 
that  the  work  described  in  thl*  report  was  undertaken. 

Our  goal  w«?  to  develop  n  capability  for  predicting  tho  fragment 
eiao  distribution  In  rock  resulting  from  dynamic  londf. .  Ke  limited  tho 
wain  of  fort  to  u  simple,  well  chornct  orlrcti,  homogeneous  rock  type, 
Arkansas  novoculi to,  under  onc-dioonaional-strain  loading.  The  approach 
consisted  of  the  following  steps: 

e  Impact  CkpcriPiontf?  wore  performed  on  u  homogeneous  rock 
under  woll-control led  conditions,  and  the  fracture  behavior 
under  various  dynamic  load  histories  wan  studied. 

e  Observations  wore  made  on  impacted  specimens  to 
establish  the  physical  processor  underlying  rock 
fragmentation. 

•  A  computational  model  of  the  f rngment otlon  process 
wss  constructed  bused  on  experimental  evidence  to 
predict  the  fragment  sizo  distribution  resulting 
from  a  given  load  history. 

•  Rock  properties  needed  for  the  model  were  measured. 

•  The  model  wn3  inserted  into  a  wave  propagation  code 
and  calculations  were  performed,'  The  rosulta  wore 
compared  with  experimental,  results. 


This  report  is  organized  Into  six  chapters.  Chapter  IX  dencrlbcie 
tho  impact  exporlniofsts,  tho  experimentol  and  calculations!  r.othcx)a  used 
to  detorElnc  the  #tros»  history  in  the  specimen©,  end  the  extent  ©f 
fracture  and  fragmentation  of  tho  BpeeiKona.  Chapter  XXI  presents  potro- 
graphic  and  fractogrcphlc  observations  made  on  the  impnesod  spwliaona; 
this  evidence  is  used  to  dodueo  the  mechanism  of  fragmentation,  i.e,,  tho 
sequence)  it  ffiieroprofce&fiea  which  occur  in  tho  tpeciaena  during  loading 
and  rosult  in  fragmentation.  In  Chapter  IV a  computational  model  of  rock 
fragmentation  under  those  loading  conditions  is  constructed  to  describe 
the  resulting  fragment  size  distribution,  Tho  model  is  based  on  the 
hypothesized  fragmentation  mechanism  and  requires  tho  loading  history 
and  certain  rock  propcrtica  as  input.  In  Chapter  V,  the  rock  properties 
that  control  tho  fragmentation  behavior  .(tho  initial  flaw  size  dlntrl- 
butloa,  tho  fracture  toughness,  and  tho  velocity  of  crack  growth)  are 
determined  in  Arkonaas  novaculite.  In  Chapter  VI  the  results  of  the 
calculations  are  presented  and  compared  with  tho  measured  recultu. 

Tho  results  of  characterization  work  on  three  other  rock  typos  ore 
presented  in  tho  ppondix. 

The  essential  rosultn  of  this  work  may  be  quickly  gresped  by 
comparing  figures  showing  experimentally  observed  crack  patterns  oed 
measured  fragment  size  distributions  for  novaculite  (figures  •!  end  S, 
respectively)  with  the  corresponding  calculntlonnl  results  shown 
respectively  In  Figures  11  ar.d  12.  Tho  ngreemont  is  considered  good, 
ond  Indicates  that  the  approach  taken  here  could  lend  to  a  useful 
solution  to  tho  problems  of  rock  fragmentation  under  dynamic  loads. 


II  IMPACT  EX  PKH I ME  NTS 


Materiel 


Arkansas  novacullte,  a  naturally  occurring  poiycryntsllina  quortrita, 
was  chosen  a a  the  bacalino  material  for  thla  atudy  because  of  its  simplicity. 
tt  is  pure,  denns,  and  Homogeneous,  and  constats  of  squlalaed,  wqulaxcd 
and  randomly  orlentod  quarts  grains  having  an  average  diameter  of  about 
10  pw.  A  population  of  flat  flaws  exists  on  planes  roughly  parallel  to 
each  other.  The  structural  and  mochenical  properties  of  novacullte  are 
presented  and  discussed  ooro  fully  in  Chapter  V. 

Impact  Loading  Procedure 

Controlled  inspect  experiments  or>  novseulito  uero  carried  out  with  a 
gas  gun  using  n  fist  projectile  impact  technique,  bo  that  fracture  and 
fragmentation  occurred  undor  one-dlmenslonal  strain  conditions. 

Projectiles,  30. 0  cm  long  by  6.33  cm  in  diameter  for  the  assallor  gas 
gun  and  76.2  cm  long  by  10.2  cm  in  diomotor  for  the  largor,  were 
accelerated  down  tho  evacuated  barrel  of  the  gun  by  the  suddon  reloose 
of  prossuro  from  nn  ndjocont  pressurized  chamber  of  helium.  Fist  tsrgot 
BpoclmoiiB  wero  impacted  with  thin  polyothyleno  or  plexiglass  flyer  plates 
attached  to  the  front  end  of  the  projectiles  (Figure  1).  Under  such  flat 

plate  impact,  comprosnlvc  wavos  initially  run  Into  the  specimen  and 

> 

projectile  head  to  produce  a  atato  of  one-dimensional  compressive  strain, 
Tonslon  Is  produced  in  thu  specimens  when  release  waves  running  Inward 
fro®  tho  free  surfneo  of  the  apoclmon  moot  similar  roloase  waves  running 
inward  from  tho  back  Hurface  of  tho  flyer  pinto.  From  a  knowledge  of  tho 
shock  Impedance  (product  of  density  and  shock  velocity)  of  the  rock  and 
flyer  plate,  the  «trer>6os  In  the  rock  apoclmon  can  be  calculated.  A 
sufficiently  large  dlnmotcr-to-thlcknoss  ratio  for  the  specimen  onsures 
that  tho  tensile  strain  is  one- diravnwlonnl  by  preventing  unloading  waves 
from  the  special*  i  periphery  from  reaching  tho  Interior  during  tho  tension 


phase. 


3 


l/*£  i  -4 


HIND 


MA-10BJ-3G 

FIGURE  1  EXPERIMENTAL  ARRANGEMENT  FOR  HIGH  RATE  TENSILE  TESTING 
OF  CYLINDRICAL  ROCK  SPECIMENS  IN  A  GAS  GUN 
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Tho  target  apecltsena  were  short  eyllndors,  usually  about  0,033  cm 
thick  nod  1,27  to  3.61  cm  In  diameter,  Their  axes  we ro  carefully  aligned 
to  coincide  with  that  of  the  gun  barrel  to  ensure  fiat  pinto  impact.  The 
speclmcna  wore  fit  tightly  Into  constraining  ringo  of  aluminum,  vhlch  in 
e  very  closo  match  to  novacullto  in  shock  impodsnee.  Most  of  the 
aluminum  constraining  rings  had  on  8~dcgrco  toper  on  the  .outer  circumfcrenco 
and  wo ro  press-fitted  Into  a  larger  aluminum  pleto  serving  as  a  epoclmon 
holder.  Others  hod  no  tapor  ond  wore  lightly  held  In  the  specimen  holder 
by  several  dabs  of  epoxy.  Upon  Impact  with  the  flyer  plate,  tho  opoclcen 
and  constraining  ring  fly  free  from  tho  aluminum  holder  one*  into  a  catcher 
tank,  which  la  filled  with  reg»  to  pi'evont  subsequent  inpaeta  and  possible 
uncontrolled  damage.  The  project llo  and  projectile  head  are  prevented 
from  entering  the  catcher  tank  by  tho  steel  plate. 

I'lf ty-throo  tensile  experiments  oh  Arkansas  novacullto  at  high  loading 
rates  were  performed  using  tho  gas  guns.  Fourteen  of  thoao  wore  instrumented 

with  ytterbium  stress  pagos  to  determine  stress  histories  or.d  to  measure 

•l 

the  magnitudes  of  rccompression  waves  produced  by  fracturo,  An  additional 

twonty  wero  uninetruraonted  experiments  in  which  wo  attempted  to  doton/iinc 

tho  dynoolc  frncture  strength  of  novncullte  and  any  effects  of  orientation 

12 

on  fracture  strength,  ’  whereas  ton  others  woro  performed  on  Jacketed 
specimens  to  study  fragmentation  behavior.  The  remaining  nine  experi¬ 
ments  wore  either  of  a  preliminary  nature,  to  determine  in  u  rough  way 
tho  impact  velocities  and  oibrr  experimental  conditions  required  to 
carry  out  tho  proposed  program,  or  of  on  investigative  nature  to  determine, 
for  example,  the  reason  for  tho  occurrence  of  undesirable  radial  crocking. 
Determination  of  Dynamic  Tenailo  Strength 

A  scries  of  20  uninstrumen  ;ed  oxporimentn  was  performed  to  measure 
tho  dynamic  tnnsilo  strength  of  novaculite,  to  determine  the  effect  of 
specimen  orientation  on  tho  dynamic  fracture  strength,  and  to  produce 
fracture  surfaces  for  examination  in  the  scanning  electron  microscope 
to  gain  information  concerning  tho  fracture  mechanism. 

3' 


Ten  specimens  wore  cut  so  ttuu  the  Impact  direction  was  normal  to  tho 
pianos  of  the  inherent  flaws;  spec incus  In  tho  remaining  10  were  oriented 
»o  that  tho  Impact  direction  was  parallel  to  the  flaw  planes.  All  speci¬ 
mens  wore  1.27  era  In  dlaraeter  by  0.035  cm  thick  and  press -fit,  tad  Into 
3.81  cm  diameter  by  0,633  cm  aluminum  onnult.  Tho  outhur  periphery  of 
each  annulus  vea*  provided  with  an  8-degreo  taper  to  fncllitute  ojoction 
upon  impact  of  tho  epoclmea-rlng  assembly  iron  the  aluminum  target  plate. 

Very  extensive  cracking,  usually  resulting  In  comminution  and  lose 
of  one  or  both  sides  of  tho  specimen  from  tho  aluminum  annulus,  was 

produced  ih  all  of  the  first  ntno  experiments  ©von  at  etreenca  a«  low  as 
2 

20  HN/n  (2000/psl).  Tills  initially  puzzling  behavior  was  prevented  In 
tho  subsequent  11  shots  by  lightly  tacking  tho  tapered  annuli  to  tho 
torgot  plates  in  throo  or  four  places  with  epoxy  Instead  'if  press-fitting. 

In  tho  lattor  oxpcrteiontP  no  radial  cracking  occurred,  and  n  good  ostimnt© 
of  t^o  dynamic  fracture  strength  wna  obtained. 

Tho  ronultn  of  this  series  of  experiments  nre  presented  In  Table  1. 

Tho  dynamic  tensile  strength  was  taken  as  tho  overage  of  the  highest  stress 
at. which  lid  damage  could  be  observed  on  diametrical  sections  of  tho 
specimen,  and  the  lowest  stress  ut  which  Incipient  apnllnt.Km  occurred. 

A  crosB-soctlonnl  view  of  Incipient  spallation  la  provided  in  Figure  2. 

For  tho  pnrallol  spocltaen  orientation,  Specimen  30  did  not.  crock  at 

’  .  2  •  2 

39, ‘I  MN/nt  (3710. psi)  but  Specimen  35  did  crock  nt  43.5  MN/m  (G300  psl); 

for  tho  normal  spocimon  orlontotion,  no  cracking  was  .observed  In  Spoclmon  43 

2 

which  was  subjected  to  n  peak  ten.clle  stress  of  38.8  MN/ra  (5030  psl), 

2 

whorens  Specimen  37  showed  cracking  nt  42.6  MN'/m  (6180  p«l).  Tims  tho 

dynamic  tensllo  strength  of  novncullte  in  the  direction  normal  to  tho  planes 

2 

of  inherent  cracks  Is  40.6  ±  2.0  MN’/ra  (5800  ±  290  psi);  the  dynamic 

2 

strength  In  tho  parallel  direction  is  41.4  i  2.0  MN/ra  (6000  ±  290  psl). 

’standard  International  (S.I.)  units,  as  now  rociuirod  by  many  government 
agencies  and  professional  soclot;ea  for  technicnl  reporln,  papers  and 
journals  will  bo  used  In  this  report.  For  convenience,  however,  English 
units  will  also  bo  given  in  parentheses  directly  thereafter.  Stress 
equivalents  to  slide  rul*1  accuracy  nre  as  follows:  1  Megapascal  (MN/'bt)  * 
1015  N'hT  a  l(! ,  2  kg/ern2,  =  145  psl  »  10  bar. 
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FIGURE  2  CROSS-SECTIONAL  VIEW  OF  SPECIMEN  30  SHOWING 
INCIPIENT  SPALLATION 
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Wo  conclude  that  tho  dynamic  tensile  tit rengt li  of  Arkansas  novneulito  io 
insensitive  to  the  orientation  of  preexisting  flows.  In  vlov  of  the 
pronounced  fin*  orientation  onisotropy  of  novneulito,  the  fracture 
strength  Isotropy  is  surprising. 

Determination  of  Fragment  Sire  Distribution 

At  low  Impact  velocities  (loss  then  cbout  20  m/soc),  Arkansas  novneulito 
specimens  could  sustain  corisidornblo  fracture  damage  while  maintaining 
integrity.  At  impact  velocities  sufficient  to  cause  fragmentation,  n 
method  was  required  to  provont  the  fragments  Oom'  scattering  and  thereby 
resulting  in  por.sible  dneutgo  and  loss  of  the  fragments,  The  arrangement 
used  for  recovering  heavily  damaged  specimens  entailed  uncasing  the 
cylindrical  rock  flpcci’vCns  completely  In  a  much  tougher  matcriul  of 
similar  shock  impedance.  Aluminum  was  found  to  he  a  suitable  encasing 
material,  first,  because  it  dobs  not •undergo  brittle  fracture  under  tho 
londli-g  conditions  of  these  experiments  and  thorafore  contains  the  cracking 
and  fragmenting  rock  specimen,  rnd, second,  because  its  shock  impedance  is 
very  similar  to  that  of  novaculito,  so  that  disturbance  of  st'reu.a  waves 
as  thoy  cross  tho  spucimon-cncasomcnt  lntcrfuco  is  minimal.'  Tho  dlmon- tons 
of  tho  turgots  were  designed  to  rcduco  edge  effects. 

As  shown  in  figure  3,  specimens  of  Arkansas  novaculito  1.27  cm  in 
diameter  by  0,6.15  cm  thick  were  fit  tightly  In  the  center  of  nn  aluminum 
disk  5.08  cm  in  diameter  and  0.953  cm  thick.  An  aluminum  cover  plato 
5.08  cm  in  diameter  by  0,318  cm  thick  was  then  placed  -over  tho  exposed 
end  of  the  specimen  and  hold  firmly  to  the  disk  with  four  equally  spaced 
screws.  An  epoxy  wna  applied  to  the  specimen  surface  adjacent  to  tho 
covor  pinto  to  ensure  Intimate  contact  with  the  aluminum  cosing.  This 
target  assembly  wao  then  subjected  to  flnt-plntc  Impact  with  tho  gas  gun. 

Ten  experiments  wore  carried  out;  the  details -are  given  In  Table  2, 

It  was  planned  to  subject  Specimens  d-J,  <15,  and  <16  to  stresses  a  factor 
of  about  2.0,  1.5,  and  1,25  in  excess  of  the  dynamic  tensile  strengths  in  un 
attempt  to  obtain  various  degreon  of  crack  coalescence  leading. to  fragmen¬ 
tation.  Tho  resulting  fracture  damage  is  described  in  tho  next  section. 
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The  next  throe  experiments,  17,  -18,  end  48  wore  performed  at  strops 
levels  near  the  dynamic  tensile  strength  to  determine  whether  the  aluminum 
encaoenont  arrangement  caused  significant  stress  amplitude  attenuation. 

If  so,  impact  velocities  sufficient  to  cause  Incipient  spall  fracture  in 
unenoasod  specimens  would  not  result  m  domago  when  encased,  Expori monte 
47  end  <19  porformed  ot  an  impact  velocity  of  10.4  m/sec  produced  signif¬ 
icant  cracking,  whereas  experiment  48  at  14.8  o/eee  produced  no  damage. 
Thcso  results  nro  in  ogreomont  with  the  damage  threshold  velocity  of 
15.1  *0,0  is/sec  eat abllmied  for  unoneased  novoaullto  in  the  first  annual 
report,*  and  so  wo  concludo  that  the  aluminum  encasement  hod  little 
attenuating  effect  on  the  stress. 

Specimens  50  and  51  were  to  bo  impacted  ot  about  23  m/scc,  In  the 
voloelty  range  of  advanced  stages  of  crack  coaloocenco  and  incipient 
t r«g.aontecion,  but  unfortunately  much  lower  velocities,  about  17  m/sec, 
wore  attained  and  much  lens  damage  resulted  than  was  dOBired.  The  final 
two  specimens  wore  shock-loaded  at  significantly  higher  velocities  to 
produce  detached  fragments.  The  aluminum  casing  (figure  3)  was  impot. .cd 
by  the  flyer  plate  in  the  flr.it  eight  oxperimcntG,  In  Experiments  52  nnd 
03,  howavor,  the  specimen  acocrably  was  oriented  nuch  that  tho  covor  ploto 
side  was  impacted. 

The  recovered  targets  were  prepared  for  f rcetographic  observation 
and  analysis  in  one  of  two  ways.  Either  they  were  cut  carefully  on  o 
dioceter  to  reveal  the  crocking  pattern  on  n  crons  section,  or  else  tho 
aluminum  oncanemont  was  removed  by  carefully  machining  the  periphery 
down  to  a  few  mils,  in  a  lathe  and  subsequently  dissolving  the  remaining 
fow  mils  of  aluminum  in  a  50«  HCL  solution. 
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In  Situ  Observations 

Specimens  45  through  52  wore  suctioned  end  polished  to  ruvool  tho 
cracking  patterns.  Tho  effect  of  stror.B  level  on  tho  extent  of  crocking 
Is  Illustrated  in  Figure  4,  which  chows  cross  ooclions  of  speciaenw 
Impacted  at  various  velocities,  Tho  choractoriatlc  dorco-shaped  crock 
pnttom  le  evldout.  Damage  Is  usually  hoavlost  In  the  half  neoror  the 
Impact  surfaco.  Fine  yarticloo  seem  to  bo  produced  at  midthickooos  and 
In  tho  roue  encompassed  by  the  dose  cracks.  L?rgo  frngmenta  originate 
mainly  near  the  flat  aurfoccs.  Tiio  froo-aurfaee  cido  of  tho  specimen  is 
usually  least  damaged  and  ie  of ton  recovered  in  one  piece,  oven  when  tho 
rcmnlndor  of  tho  specimen  has  fragmented.  Freo,  tuicoalosced  crack  tips 
nro  commonly  observed  In  specimens  Impacted  at  high  ob  well  ob  at  low 
stre&aoG. 

Fragment  Size  Distribution 

Wo  attempted  to  determine  the  frngnont  size  distributions  produced 
In  Experiments  44,  45,  and  53  by  carefully  removing  the  aluminum  onenoo- 
monto,  Specimen  44,  however,  fell  apart  In  only  a  few  largo  plecen  and 
was  unsuitable  for  n  sieve  analysis.  Specimen  45  remained  intact  after 
removal  of  the  aluminum  and  retained  consldoroblo  strength  (firm  hard 
pressure  was  insufficient  to  break  it  up),  so  it  wuo  mounted  in  epoxy 
ond  sectioned  as  described  In  tho  previous  discussion. 

Tho  fragment  size  distribution  for  Specimen  53  was  determined  by 
plocing  the  colloctod  fragments  in  tho  top  sieve  of  a  series  of  U.S. 
sieves  placed  in  tho  following  order  from  top  to  bottom:  No.  6,  tf,  10, 

14,  20,  40,  50,  100,  200,  and  400,  and  8  pan  to  enteh  tho  finos.  The 
system  wnc  vibrated  for  n  short  time,  and  the  particles  retained  on 
each  screen  were  counted  and  weighed.  Figure  5  shows  the  shapes  of  the 
particles;  the  raw  data  are  presented  in  Table  3  ond  in  Figure  G,  Then© 
experimental  fragmentation  (Intnseri  used  to  develop  and  check  tho  ,’yncmic 
froctui-o  modol. 
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FIGURE  4  POLISHED  CROSS  SECTIONS  OF  ARKANSAS  NOVACULITE 
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FIGURE  5  PHOTOMICROGRAPHS  OF  VARIOUS  SIZED  FRAGMENTS 
FROM  EXPERIMENT  S3 


Stress  Wove  Measurements 


la  fourteen  experiments,  on  attempt  woo  made  to  measure  otreoa  wave 
profiles  transmitted  through  tho  specimen. ^  in  addition  to  providing  the 
ponk  stress  and  stross  duration  experienced  by  the  rock  under  Impact  loads, 
measured  stress  histories  yield  important  information  about  tho  constitutive 
relations  and  rate  dependence.  In  the  caso  of  n  specimen  undergoing  frncturo, 
recomproasioit  wavos  emitted  as  micrbcracks  form,  lmpingo  on,  end  reload 
the  gaga  to  an  extont  proportional  to  the  dynamic  tensile  etrength  of  tho 
rock.  Furthermore,  the  slopo  of  the  reloading  pulse  gives  an  indication 
of  tho  rate  of  fracture—n  sharp  rise  corresponds  to  brittle  behavior. 

Vtterblum  piezoroaietant  stress  gages  mounted  in  plexiglass  blocks 
wore  hold  in  contact  with  the  rear  surface  of  the  specimen.  Upon  impact, 
compressive  waves  impinge  on  the  gage  and  the  stress-induced  electrical, 
signals  nro  recorded  with  an  oscilloscope.  The  tensile  stresses  In  tho 
rock  specimen  arc  then  calculable  from  tho  gage  record,  if  tho  relative 
shock  impedance  of  specimen  and  backing  taaterlal  are  known. 

The  actual  oscilloacopo  record .  obtnined  in  theso  experiments ,  as 
well  as  a  more  detailed  account  of  this  stross  wave-measuring  technique, 
can  be  found  in  reference  1.  The  results  may  be  summarized  by  noting 
that  clear  and  consistent  differences  in  the  gage  records  were  obsorvort 
for  specimens  which  underwent  fracture.  Therefore  it  was  possible  to  toll 
whether  or  not  a  specimen  was  cracked  before  actually  examining  it.  Peak 
compressive  stresses  experienced  by  the  specimens,  however,  were  not 
mensurable  from  the  gage  records  because  of  significant  deviations  from 
planar  impact  and  hence  excessive  rounding  of  the  oscilloscope  traces. 
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Ill  .THE  MECHANISM  OF  FRAGMENTATION 

Considerable  importance  la  attached  to  establishing  tho  fragmentation 
mechanism,  because  this  Tonus  tho  basin  Tor  the  computational  modal.  The 
fragmentation  mechanism  in  defined  an  that  sequence  of  event*  that  occurs 
in  the  rock  during  loading  nnd  leads  up  to  f rogmontation.  We  constructed 
tho  computational  model  for  predicting  fragment  size  distributions  by 
modoiing  individually  the  physical  processes  that  precede  fragmentation. 

A  model  developed  on  such  a  physical  basis  should  bo  inherently  ,aore 
accurate  an  roll  as  more  applicable  to  other  materials  under  other  loading; 
conditions  thn^i  model  based  on,  nay,  empirical  correlations.  Thus,  eerh 
impactod  rock  spocimon  was  exnmlnod  carefully  in  search  of  evidence  con- 
corning  tho  physical  processes  involved  in  fragmentation. 

Froctographlc  Techniques 

Two  main  f ractographic  techniques  woro  used  to  deduce  the  fragmentation 
mechanism  from  impacted  specimens.  In  one  technique  fracture  surfaces  nnd 
ln-’.vidual  fragments  wero  examined  by  optical  and  scanning  electron  micro¬ 
scopy  to  look  for  murkings  that  would  reveal  how  the  surfaces  nnd  fragments 
woro  formed.  Those  observations  providod  information  on  tho  histories  of 
Individual  cracks.  The  other  technique  entailed  Btopping  the  fragmentation 
process  at  various  stages  of  completion  through  careful  control  of  the 
experiment,  then  sectioning  the  specimen,  and  oxnmining  tho  pattern  of 
cracks  intersecting  tho  section.  This  technique  yieldod  insight  into 
the  behavior  of  interacting  cracks. 

Fractogrnphic  Observations 

The  purpose  of  optical  nnd  scanning  electron  microscope  examinations 
of  fracture  surfaces  was  to  look  for  inhomogeneities  that  could  have  served 
os  crack  initiation  sites,  and  river  lines  and  hesitation  lines  which 
indicate  the  nature  of  crack  growth. 


As  la  typical  of  most  rocks,  tho  fractured  surfaces  were  nearly 
featureless  and  for  the  most  part  yielded  llttlo  evidence  of  how  fracture 
occurred.  The  broken  surfaces  consisted  of  countless,  well-defined  and 
o^ulslzed  polygonal  blocks— quartz  groins  exposed  by  tho  passage  of  > 
crack  along  the  grain  boundaries,  Only  Infrequently  were  propagation 
markings  and  heal tat Ion  lino*  observed.  The  circular  markings  In 
Figure*  7b  and  b  ere  know  as  hesitation  lines  because  thoy  are  produced 
when  the  crock  undergoes  *  sudden  change  of  velocity  as,  for  Instance, 

"  wheA  ! the  crack  id  impinged  upon  by  «  utresa  rev*.  They  are  analogous  to 
the  orreat  lines- formed  where  a  crack  he»  actually  stopped.  Thus 
hesttotion  lines  and  arrest  lines  delineate  the  position  of  tho  crack 
front  at  some  instant  in  time  and  thereby  reveal  tho  contour  of  the 
propagating  crack.  From  the  .appearance  of  the  hesitation  lines  in 
Figure  7,  wo  deduce  that  the  cracks  had  circular  peripheries,  i.e.,  they 
were  penny-shaped, 

Alsooviderit  on  the  fracture  surface  in  Figure  7a  arc  several  lines 
radiating  outward  from  the  apparent  center  of  the  circular  hesitatio" 
lino  and  intersecting  the  hesitation  lino  ot  right  angles.  These  lines 
are  actually  height  discontinuities  and  thus  are  rominlscent  of  river 
lines  and  cleavage  steps  commonly  observed  on  fracture  surfaces  of 
metallic  and  lonlcully  bonded  materials.  As  such  they  form  parallel  to 
the  propagation  direction  and  thus  indicate  the  propagation  direction. 
Thus  the  river  lines  provide  a  running  history  of  the  crock  path  and  n 
hesitation  lino  gives  the  crock  shape  at  an  instant  in  time.  Figure  7 
Is  therefore  interpreted  to  mean  that  crack  proportion  in  novacvllto 
under  dynamic  uniaxial  strain  loading  occurs  by  tho  expansion  of  penny¬ 
shaped  cracks  radially  outward  from  initiation  sitos. 

A  corollary  of  the  fact  that  river  lines  form  parallel  to  tho  crack 
propagation  direction  is  that  the  fracture  origin  may  be  located  by 
following  the  linos  in  a  direction  opposite  to  the  propagation  direction. 


Figure  7  shows  clearly  that  the  river  linen  emanate  from  n  slnglo  site 
situated  at  the  approximate  center  of  the  circular  hesitation  line. 

High  resolution  examination  of  this  site  with  the*  scanning  electron 
microscope  did  not  reveal  any  obvious  heterogeneity  that  might  have 
servod  ns  the  weak  spot.  In  fact  in  no  instance  was  it  possible  to 
identify  positively  the  hotorogenoity  in  the  rock  responsible  for  the 

initiation  of  a  crack  at  that  particular  place.  Thus  wo  can  only 

0 

speculato. that  the  heterogeneities  responsible  for  crack  Initiation 
were ‘the  small  crack-like  flaws  that  existed  inherently  In  the  rock. 

Tho  nature  of  these  flaws  is  described  mo.ro  fully  in  Chapter  V. 

Individual  fragments  such  as  those  shown  in  Figure  3  woro  also 
examined  with  optical  microscopes.  Their  shapes  ns  con  be  seen  from 
tho  figure  arc  roughly  ccjulaxed,  as  opposed  to  olongnted  bodies  us 
might  perhaps  be  expected  from  Figures  4b  nnd  4c.  In  general  tho  fragments 
had  from  6  to  8  facets,  regardless  of  fragment  size,  implying  that  on 
the  average  3  to  <1  cracks  are  associated  with  ono  fragment,  and  indicating 
that  larger  cracks  produce  larger  fragments  whereas  smaller  cracks 
prod'  .  a  smaller  fragments.  It  was  also  noted  that  individual  fragments 
generally  contained  a  number  of  cracks  in  their  volumes,  which  means 
that  not  nil  cracks  were  effective  in  producing  fragments. 

The  purpose  of  performing  impact  experiments  at  stresses  insufficient 
to. produce  fragments,  was  to  provide  specimens  for  observation  that 
contained  various  degrees  of  prof rngncntntlon  fracture  damage.  Such 
specimens  woro  cut  carefully  on  n  diameter  with  n  diamond  saw,  and  this 
surface  was  polished  to  reveal  n  cross  soction  of  the  crocicing  pattern 
and  to  provide  in  situ  views  of  fragmentation.  Figure  4  shows  the 
nppournneo  of  four  such  specimens  impacted  at  stress  levels  ranging  from 
one  below  the  tensile  strength  through  one  sufficient  to  produce  loose 
fragments.  This  figure  may  also  be  viewed  as  depicting  successive  stages 
in  time  of  a  fragmenting  rock  specimen,  i.e,,  ns  the*  stress  in  the  wave 
rises,  increasingly  more  cracks  originate  and  grow. 
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Those  pictures  show  that  crocks  Initially  tom  nonr  to  the  impact 
surface  and  propagate  predominately  in  the  lateral  direction  (normal  to 
the  direction  of  tho  maximum  tonsiic  stress).  At  low  stresses,  lew 
fragments  can  torn,  and  thoy  arc  large.  As  the  stress  level  increases, 
more  crocks  ure  initiated  and  moro  coalescence  occurs.  Crock  growth  does 
i •  >t  appear  to  be  particularly  sensitive  to  stress  level.  Only  at  tho 
highest  stress  (Figure  4d)  can  significant  numbers  of  vortical  cracks  be 
seen;  Tho  effect  of  the  vertical  cracks  is  clearly  to  increase  the  extent 
of  specimen  cosuainution.  They  transfom  the  large  elongated  fragments 
(Figure  4c)  to  many  smaller,  more  equloxed  fragments  (Figuro  4d)  and  thus 
influence  strongly  the  resulting  frag;  .nt  size  distribution. 

Figure  4d  also  indicates  that  smaller  fragments  tend  to  originate 
near  the  raid-thickness  of  the  specimens  whorons  the  larger  fragments  came 
from  the  motorlnl  near  tho  impacted  and  free  eurfnr  »s.  . 

t 

Hypothesized  Fragmentation  Mechanism  >  . 

.  From  those  observations  wc  envision  a  mechanism  '<  i  fragmentation 
consisting  of  four  stages,  namely: 

(1)  Activnt  on  of  a  number  of  preexisting  structural  flaws 

(2)  Propagation  of  activated  cracks  radially  outward 
■  (3)  Coalesccnco  and  branching  of  propagating  crocks 

(•!)  Isolation  of  individual  rock  fragments  from  one  another. 

Tho  first  two  st.ngos  have  been  verified  by  similar  impact  experiments 
on  a  transparent  material.  Figure  8  shows  the  internal  penny-shupod  cracks 
produced  in  this  material,  a  polycarbonate.  The  tiny  blnck  areas  that  are 
obser’".ole  at  tho  center  of  every  crack  have  been  positively  identified  as 
Inherent  flaws  in  the  material  produced  during  fabrication.  The  multiple 
concentric  rings  visible  oil  the  crack  surfaces  in  polyearbonato  arc  similar 
to  the  Hesitation  lines  found  on  novaculitc  fracture  surfaces  and  imply 
that  the  crack  morphology  in  novaculitc  should  be  basically  penny-shaped 
and  should  remain  penny-shaped  during  growth  until  coalescence-  begins. 
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TV  CONSTRUCT  ION  01'  TS1E  PTlEDlCTIVK  CAPABILITY 


Tho  capability  for  predicting  the  fragment  size  distribution  for 

novaculite  under  dynamic  ono-dimensioriEil-straln  loads  was  obtained  by 

combining  a  computational  model  ot  the  fragmentation  mechanism  with  a 

finite  difference  wave  propagation  code.  Tho  computational  fragmentation 

model  is  based  on  tho  f ragnontatiou  mechanism  hypothesized  in  tho  provlouo 

chapter,  and  treats  quantitatively  each  of  the  four  stages,  namely; 

flaw  activation,  crack  growth,  creek  coalescence^  and  f t*if;Hent  formation, 

4 

The  one-dimensional  wave  propagation  code,  PUFF,  calculates  from  tho 
impact  conditions  the  stress  history  in  the  specimen,  l.e, ,  tho  magni¬ 
tude  and  duration  of  the  stress  at  any  point  in  time. 

Tho  Computational  Fragmentation  Model 

The  fragmentation  of  novaculite  under  dynamic,  onc-dtmens tonal- 
strain  conditions  is  modeled  by  treating  quantitatively  tho  hypothesized 
stages  of  tho  fragmentation  mechanism,  A  treatment  of  each  of  the  four 
et  -jos  la  given  below. 

Activation  of  prooxlstlng  flaws — Based  on  experimental  observations 

of  impacted  and  tin  impacted  Arkansas  novacullto  (Chuptors  III  and  V, 

respectively)  we  assume  that  the  matorial  contains  inherently  a  population 

of  pcnny-slmped  flaws  of  varying  sizes,  some  of  which  become  unstable 

upon  passage  of  the  stress  wave  end  develop  into  propagating  cracks.  To  . 

calculate  the  number  of  flaws  activated  by  a  given  stress  pulse  (and 

hence  tho  number  of  cracks  in  the  material)  wo  invoke  r»  Griff ith-Irwin 

5 ,  G 

fracture  mechanics  criterion,  i.o,,  flaws  having  a  radius  c  larger 

than  some  critical  size  that  is  a  function  of  the  tensile  stress 
*  ♦  ♦ 
c  =>  c  (cr)  will  be  activated,  whereas  those  flaws  smaller  than  c  remain 

dormant.  Of  tho  various  crack  geometries  for  which  fracture  mechanics 

7 

stress  analyses  exist,  Sneddon's  relation  for  an  internal  penny-shaped 
crack  in  an  infinite  elastic  medium  subjected  to  n  uniform  tensile  stress 
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normal  to  tho  crack  piano  i 

•  22  * 

C  =  1TK  /4TT  (1) 

Ic  •  : 

i 

rnoBt  nearly  applies  to  tho  conditions  of  our  experiments,  i.e.,  tho  5 

cracks  appear  to  bo  roughly  penny-shaped,  and  tho  specimen  behavoB  ns  • 

•  .  *•  < 
if  it  wero  infinite  during  the  life  of  tho  tensile  pulse.  Tho  piano  * 

strain  fracture  toughness  K,  for  novaculite  ia  determined  In  tho  !’ •• 

Ic  i  . 

following  chapter.  Thus  wo  uso  Kq,  (1)  to  calculate  the  size  of  tho  IY. 

..  ‘  *  .  '  .  ‘  •  .  .  '  1-  ;  ’ 

amallest  flaw  which  will  bo  activated  by  a  dynamic  load.  To  obtain  the 

•  ' 

number  of  activated  flaws,  wo  need  to  know  tho  slzo  distribution  of  flaws  > 

In  tho  rock.  This  has  been  measured  for  novaculite  and  tho  results  are 

i 

presented  In  Chapter  V. 

Crack  Propagation — The  distance  which  each  crack  can  propagate 

depends  on  the  crack  velocity  and  tho  duration  of  the  stress  pulse.  Wo  : 

assume  here  that  the  cracks  accelerate  very  rapidly  to  a  constant  maximum 

velocity  of  one  third  the  longitudinal  wave  speed.  This  Is  in  accord 

with  our  experimental  observations  (Chapter  V)  cs  well  as  with  theoretical 
1 1 

estimates.  Knowing  that  crack  propagation  occurs  radially  outward 

from  activation  sites,  tho  fracture  surfaco  area  produced  per  crack  nt  ; 

uny  given  time  step  can  be  approximated  by 

7 

A  (t>  a*  2ff(C.&t/3) 
i  1 

providing  tho  cruck  was  not  stopped  prematurely  by  barrier  to  crack 
growth  or  coalescence  with  other  cracks.  Individual  cracks  arrest  when 
the  stress  level  falls  below  the  value  given  by  Kq.  (1).  Thus  the 
total  fracture  surface  areu  produced  in  novaculite  by  a  known  stress  pulse 
of  peak  stress  (7  and  duration  £t  insufficient  to  cause  significant 
coaloscenco  Is  given  by 


(2) 
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total  J  l 


where  N  is  tho  number  of  Activated  flows, 
net 

Furthermore,  tho  total  energy  absorbed  in  creating  new  surface 
area  can  bo  estimated  If  the  specific  fracture  surface  energy  V  i» 
known. 


E#  '  .  "  v 

fracture 
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Crack  Coalescence— Our  quantitative  treatment  of  coalescence  of 
p.-opugating  cracks  ontailod  assuming  a  criterion  for  tho  distance 
between  converging  cracks  at  which  interaction  first  occurs,  and  taking 
a  statistical  approach  by  treating  largo  numbers  of  coalescing  cracks. 

To  develop  tho  coalcscenco  criterion,  tho  concept  of  the  "crack 

rango”  is  introduced.  Tho  crack  range  T  (c)  is  dofined  as  the  volume 

c 

of  solid  material  surrounding  a  crack  that  experiences  magnified  strains 
of  some  arbitrary  level,  which  arise  from  the  stress  concentrating  effect 
of  the  crack.  When  tho  crack  ranges  of  two  cracks  overlap,  they  are 
considered  to  have  coalesced.  As  an  upper  limit (obtained  from  the 
solution  of  tho  elastic  stress  field  around  a  penny-shaped  crack  in  an 
infinite  modiura)  tho  crack  range  includes  the  material  within  one  crack 
radius  of  uny  part  of  the  crack.  Thus  for  a  penny-shaped  crack  tho 
maximum  crack  range  Is  an  ellipsoidal  volume  of  revolution  shaped  like 
a  hamburger.  The  volume  then  Is 


I 


(  .» 


i’  "V, 
) 

•  :4 

i  • 

i 

>  ■; 

i 

■  . 
I  t 


\ 

i  ■  • 


2  TO  3  3 

T  (c)  a  (IT  +  —  IT) c  ~  20.34  e 
C  3 


(5) 


The  observation  that  the  fragments  are  nearly  equlaxed  suggosts  u 
much  smaller  crack  rnr.ge  volume.  For  lnstanco  if  the  fragments  were 
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cubo-»h.ipod ,  then  the  crack  ranges  would  bo  tho  volumes  of  two  small 
pyramids  on  either  aide  of  onch  crack  (ono  crack  forming  ono  ride  of 
tho  cube)  and  tho  crack  range  volume  is 

,  ;  •  r  (o)  «  |  c3  (6) 


For  our  purposes  wo  will  define 


T 

c 


by 


T  (c)  «  T  o3 

c  -  ■;  ,  c  .•> 


where  T  will  probably  lie  between  2  and  20. 
c 

We  con  now  poatulato  a  coalescence  criterion.  It  Is  logical  that 
two  propaguting  cracks  in  an  infinite  body  hRvo  no  knowledge  of  one 
another  until  they  come  within  some  critical  distance  of  each  other. 

When  this  critical  nearness  ia  attained,  the  cracks  sense  each  other 
and  interaction  begins.  The  concentrated  stresses  at  each  crack  front 
superimpose,  and  the  propagation  behavior  of  oBCh  crack  is  influencod 
by  tho  presence  of  the  othor.  In  tho  present  experiments  the  crack  tip 
stresses  are  similar  and  therefore  should  be  additive,  thereby  encouraging 
the  cracks  to  propagate  into  each  other  and  coalesce.  The  extent  of  the 
strain  field  about  a  crack  is  given  by  its  rnngo  as  previously  defined, 
and  wo  proBuino  that  crack  coalesccnco  occurs  when  the  ranges  overlap. 

For  numerical  computations  with  tho  wave  propagation  code  we  adopt  tho 
following  coalosconco  criterion:  crack  coalescence  occurs  when  tho  sum 
of  tho  crock  rango  volumes  for  all  tho  cracks  in  a  finito  difforence 
cell  is  equal  to  tho  volume  of  that  coll.  This  criterion  may  be 
expressed  mathematically  as  follows: 


Consider  a  distribution  of  penny-shaped  cracks  whose  size  distribu¬ 
tion  is  given  by  p  (c),  tho  density  oi  cracks  as  a  function  of  size, 
c 
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Here  wo  noto  that  the  fragment  size  distribution  is  a  Junction  of  tho 
radii  of  tho  cracks,  not  of  tho  fragments.  Lot  Tf(cf)  bo  t:.o  volume 
of  tho  fragment  with  fragment  radius  cf.  This  volume  nay  be  written 


T  (c  )  T  c 
f  f  f  f 


(12) 


whoro  T  is,  for  example  4n/3  for  spherical  particles  and  8.0  for 
cubes. 

The  fragment  volume  will  just  fill  tho  total  volume  at  the  t<mo  of 
complete  fragmentation.  Therefore  tho  fragmentation  critorlon  is 
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To  determine  tho  number  and  size  of  fragments,  wo  introduce  the 
following  relations  between  crack  nnd  fragment  variables. 


o'  «  v  Bfl  (c) 
s  r  c 


(13) 


(14) 


(15) 


cf  “  Y  c 


(16) 


where  8  nnd  y  may  be  functions  of  c,  and  V  is  the  remaining 
relative  volume  which  has  not  been  fragmented  by  cracks  larger  than  c. 
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For  8-sidod  fragments  whore  each  face  is  formed  by  one  crack.  S  a  1/4 
and  y  =  1.0, 

The  procedure  for  computing  tho  fragment  Blase  distribution  and 
for  determining  whether  fragmentation  has  occurred  begin*  with  V  «  1.0. 

T 

First,  tho  largest  cracks  arc  transformed  to  fragment#  using  Eq.  (15). 
These  fragments  have  a  volume 
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da  t  y  c 

i  f  i 


(17) 


in  tho  computer  program.  This  volume  V  is  presumed  to  contain 
small  cracks  which  do  not  form  fragments.  Therefore,  to  compute  tho 
fragments  for  the  next  smaller  crack  sizo,  thin  volume  is  removed  bofora 
computing  tho  number  of  fragments  of  tho  next  smallor  size,  ■ in  general, 
the  lumbor  of  fragments  is 
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l  f 

where  •  V  «  1  *-  J  y  , 
r  J=U1  J 

This. process  for  determining  the  number  of  fragments  is  followed  until 
nil  activated  cracks  nru  treated  or  until  nil  the  volutno  is  fuli  cf 
fragments.  If  nil  the  volume  is  filled  with  fragments,  then  tho  cell 
is  presumed  to  be  completely  fragmented. 

At  tho  end  of  tho  computation  the  material  that  contains  coalescing 
cracks  but  is  not.  fully  fragmented  is  taken  as  partially  fragmented. 

The  number  of  fragments  of  each  size  Is  determined  by  the  foregoing 
calculation.  The  remaining  unfrngmented  material  is  not  assigned  to  any 
fragmont  size  group,  but  is  assumed  to  bo  part  of  one  of  the  large  chunks 
observed  in  tho  experiments . 
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S r rosa  History  Calculations 

4 

FUFF ,  #  one-dimensional,  finite  difference  wave  propagation  code 
was  used  to  calculate  the  Magnitude  and  duration  of  the  stress  at  ovory 
location  in  the  spociuten  and  at  every  point  in  tirao.  It  also  calculates 
tho  current  density  and  energy.  Because  of  tho  low  stresses  required  to 
produce  fragmentation,  the  Kugoniot  elastic  limit  is  never  exceeded  and 
plastic  behavior  Is  never  realised.  Therefore  a  very  slmplo  equation 
of  state,  Hooke’s  law,  can  be. used  to  calculate  the  st.roB603  in 

12 

novaculito.  The  equation  of  state  for  plexiglass  is  well  known.  The 
PUFF  code  has  bocm  used  successfully  to  predict  wave  profiles  in  a 
number  of  materials  under  known  conditions  of.  one-dimensional  impact 
(Reference  13,  for  cxamplo). 

The  simple  elastic  equation  of  state  for  novaculito  is  no  longer 
adequate  for  describing  the  material  response  once  the  fragmentation 
process  begins.  When  cracks  begin  to  form  and  grow,  the  apparent 
elastic  stiffness  decreases,  fiocompression  waves  run  out  from  the 
surfaces  of  the  propagating  croc  :  to  Interact  with  and  orode  the  stress 
pulse.  This  effect  is  of  considerable  magnitude,  and  must  bo  taken  into 
account  in  the  calculations.  ■  ’ 

Tho  procedure  used  to  dctornlnc  tho  stiffness  of  material  under- 

13 

going  fracture  and  fragmentation  has  been  given  elsewhere.  The 

basis  of  the  method  is  the  concept  of  a  two-component  system:  solid 

material,  and  void  inside  the  open  cracks.  Tho  specific  volume  of  the 

system  changes  when  loaded:  this  change  is  due  partially  to  tho  solid 

and  partially  to  the  chunge  in  the  void  volume.  The  void  volume  calcula- 

7 

tion  follows  the  analysis  of  Sneddon  for  a  penny-shaped  crack  in  an 
elastic  material.  During  the  wave  propagation  calculation,  the  following 
sequence  of  events  may  occur  while  the  material  is  in  tension: 
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(a)  Initial  tensile  loads  cause  only  olastlc  volume  changes 
in  the  solid  until  u  threshold  stress  (equal  to  the 
dynamic  fracture  strength)  is  reached. 

(b)  Whon  the  stress  exceeds  the  frueture  threshold,  cracks 
begin  to  nucleate  £u;d  grow,  and  tho  void  volume  produced 
by  the  cracks  acts  to  decrease  tho  volume  change  required 
o£  tho  solid.  Thus  tho  tensile  stress  in  tho  solid  is 
lower  than  it  would  bo  .or  undamaged  material  under  tho 
same  volumo  change;  henco,  the  effective  modulus  of  tho 
solid  has  decreased. 

(e)  With  increasing  stress  and  continuing  volume  change,  a 
point  is  reached  whore  the  void  volume  increase  just 
equals  tho  applied  volume  change.  Hero  there  is  no 
change  in  solid  volume  and  henco  no  change  in  tensile 
stress:  tho  stress-volume  path  has  ruoched  a  peak  and 
tho  effective  modulus  is  zero. 

(d)  With  further  volume  changes  the  increase  of  void  volume 
(by  growth  and  nucleation  of  cracks)  tends  to  exceed  tho 
applied  volume  change.  Then  the  solid  volume  change  is 
negative  and  the  tensile  stress  is  decreasing.  The 
effective  modulus  is  negative  during  this  period. 

Co)  If  compressive  volume  changes  occur  at  any  time,  there 
will  usually  bo  a  decrease  in  both  void  volume  and  solid 
volume.  Then  the  effective  modulus  is  positive  as  in  (a). 

(f)  If  no  net  volume  change  occurs,  there  will  usually 
be  an  increase  in  void  volume  and  a  corresponding 
decrease  in  solid  volume:  the  effective  modulus 
in  such  a  case  is  infinite. 

The  preceding  sequonce  of  events  is  treated  in  the  proaont  formula¬ 
tion  of  tho  behavior  of  materiel  undergoing  fracture,  A  derivation  of 
tho  effect  ol’  damage  on  the  material  stlffnoss  is  glvon  below  to  indicate 
the  basis  of  tho  method.  Hero  the  case  is  considered  in  which  no 
growth  or  nucleation  of  cracks  occurs:  the  resulting  modulus  is  that 
which  is  appropriate  for  a  residual  strength  calculation.  The  funda¬ 
mental  relation  is  that  the  total  volumo  change  is  the  sum  of  solid  and 
void  volume  changes: 
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£V  =  LV  +  &V 
a  v 


( 1.9) 


7 

The  void  volume  change  ia  dorlved  from  the  analysis  of  Sneddon  for  the 
opening  of  «  ponny-shapod  crack  in  an  elastic  medium  under  uniform 
tension.  The  half-opening  of  tho  crack  faces  is 

„  4(i-v2)ccr 

o  «* - - - 

rrK 

•  where  c  is  the  half  crack  length 

a  in  tho  stress  applied  normal  to  the  crack  plane 
E  is  Young’s  modulus 
v  is  Poisson’s  ratio. 

The  void  within  the  crack  faces  is  an  ellipsoid  with  semlaxos  c,  c,  and 
6,  so  tho  volume  is 


V 

c 


16(l-v2)  c3(J 


(20) 


To  determine  tho  total  crack  volume  n  sum  is  made  over  all  cracks. 


whore  N  is  tho  number  of  cracks  of  volume  V  per  unit  volume.  Then 
1  ic 

Eq.  (19)  becomes 


(21) 
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where  Let  U  a  change  in  stress 

K,  u  are  bulk  and  shear  moduli 
V  ia  the  total  specific  volumo. 

Dividing  Eq.  (21)  by  £o ;  an  ciffective  modulus  M  is  dorlved. 
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K  +  -U 

This  compliance  ie  similar  td  that  which  would  bo  obtained  for  a 

4 

composite  made  of  solid  material  plus  void.  Tho  term  K  +  7  p  ia  the 

J 

stiffness  of  the  solid  and  V  /V  is  its  volume  fraction,  Tho  last  term 

s  . 

on  the  right  in  Eq,  (22)  must  be  into.-pretod  then  as  the  compliance  of 
the  cracks.  This  expression  shows  that  the  compliance  of  the  composite 
increases  with  increasing  void  fraction,  l.e. ,  Increasing  number  and/or 
Increasing  size  of  cracks. 


*  '  i 

I  i 


i.  1 


i  * 

v  .  • 


35 


•tfl'.t'iw OMlAt  y»*w  ***• Wk.. 


V  HOCK  PKOPKltTIES 

ti<ccutton  of  the  predictive  capability  described  in  the  provious 
chapter  requires  a  knowledge  of  certain  rock-apecif ic  properties,  for 
their  values  are  required  as  input.  Included  are  the  size  distribution 
of  inherent  flaws,-  tho  plane  strain  fracture  toughness,  and  the  crack 
growth  velocity.  Values  of  theso  parameters  for  Arkansas  novaculite  are 
determined  in  tjis  chapter. 

Inheront  Flaw  Structure  -  •*  .  ■  ...  .. 

m  n~  **"*  '  r  1  1  "■  V *■  *  r  H. f.  if  ,  -  ‘Y  •■iju”*'.  *!•  '**•''  *’ 

•  '  '  ,  ■  •  .  • *-4*:  V  ;  -  V:  •' 

Casual  observation  with  the  unaided  eye  of  polished  novaculite  surfaces 
Is  sufficient  to  show  thit  this  rock  has  a  strongly  oriented  flaw  structure. 
The  flaw  structure  is  clearly  evident  in  the  low  magnification  composite 
micrograph  of  Figure  9.  .The  size  distribution  of  tho  Inherent  flaws  was 
determined  by  counting  and  measuring  the  flaw  traces  on  a  polished  section 
through  a  specimen.  Those  data,  which  represented  the  size  distribution 
of  flaw  traces  per  unit  area,  were  then  converted  by  means  of  a  statistical 
transformation  to  obtain  the  actual  size  distribution  of  inheront  flaws 
po.  unit  volume. 

Nino  overlapping  photographs  at  40X  wore  required  to  span  the 

diameter  of  n  rock  specimen.  Tho  total  photographed  area  was  slightly 
2 

more  than  1.1  cm  ,  and  in  this  area  19-1  preexisting  crnck  traces  were 

counted  and  measured.  These  data  wore  converted  by  means  of  a  Scholl 

14  15 

typo  statistical  transformation  implemented  by  tho  BABS2  computor 

16 

code;  tho  results  are  presented  in  Figure  10. 

More  tho  cumulative  concentration  of  cracks  having  radii  greater 

than  radius  c  is  plotted  ns  a  function  of  c  .  The  relatively  few  (194) 

1  2  * 

traces  observed  on  about  1  cm  of  the  surface  of  section  transform  into 

3 

a  very  large  volume  density  (~1 00 , 000/cm  ).  The  size  distribution  of 
preexisting  crack3  hus  a  parabolic  form  in  log-normal  space  with  a  cut¬ 
off  at  about  c  =  500  urn.  The  curve  is  well  described  by  the  analytical 
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FIGURE  9  COMPOSITE  MICROGRAPH  OF  A  BLOCK  OF 
NOV ACU LITE  SHOWING  THE  PREFERRED 
ORIENTATION  OF  THE  INHERENT  FLAWS 


expression 

N  =  6(c)  exp[ll . 1  -  (3.7  x  102)c  +  (0.42  x  104)C2]  (23) 

whero  6(c)  a  l  for  0  <  c  £  0.03 
6(c)  =  0  for  c  >  0.05 

and  c  is  in  centiinoterw. 

The  cutoff  in  crack  half  siz.0  at  about  500  um  appears  to  be 
realistic.  A  very  large  number  of  observations  of  polished  novaculitc 
surfaces  wore  made,  and  nover  were  flaw  traces  significantly  greater 
than  1  mm  observed.  Statistically  speaking  if  many  observations  of 
crack  traces  ore  made,  the  length  of  tho  largest  trace  is  approximately 
equal  to  the  diameter  of  the  largest  flaw.  Thus  the  radius  of  the  largest 
inherent  flaw  in  the  present  specimens  of  Arkansns  novaculito  was  taken 
to  be  about  500  pm. 

Fracture  Toughness 


The  plane  strain  fracture  toughness  K  is  a  material  property 

Ic 

that  describes  the  resistance  of'  tho  material  to  crack  propagat J on ,  This 
parameter  is  used  in  the  first  stage  of  the  fragmentation  model  to 
determine  the  critical  flaw  size  for  a  given  dynamic  stress  and  hence 
the  number  of  inherent  flaws  in  novaculito  that  became  propagating 
cracks.  . 

Tho  dynamic  plane  strain  fracture  toughness  for  novaculito  was 

2 

calculated  from  the  measured  dynamic  tensile  strength  of  41.0±2.0  MN/m 

(5950  ±290  psi)  (Chapter  II)  and  the  radius  of  the  largest  flaw,  500  um, 

established  in  tho  previous  section,  Sneddon's  expression,  Eq.  (1),  is 

usod  to  relni e  tho  dynamic  tensile  strength  and  the  radius  of  the 

largest  penny-shaped  flaw  c  to  the  dynamic  fracturo  toughness  K 

max  Ics 


i 
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and  a  value  for  K  of  1.04  MN’/rc  is  obtained.  Quasi-static  tensile 
Ics  j 

strength  determinations  using  the  SR!  expanded  ring  tests  yielded 

similar  values  for  the  tensile  strength  and  hence  the  fracture  toughness, 

demonstrating  that  these  properties  of  novaculite  are  strain  rnte- 

— i  -l 

insensitive,  at  least  for  strain  rates  in  the  range  of  .3  x  10  sec 
3  '--a 

to  7  x.  10  sec  . .  >».  - 

Crack  Velocity 

The  distance  each  crack  can  propagate  depends  on  the  crack  vulocitj- 

and  the  duration  of  the  stress  pulse.  Ytterbium  stress  gages  wero  used 

to  measure  the  latter  in  this  work  and  crack  velocities  were  inferred 

from  direct  measurements  of  the  distance  propagated  by  cracks  in 

essentially  crack-free  material.  The  radii  of  the  cracks  shown  in 

Flguro  2  are  about  2  mm.  Tho  ytterbium  stress  gftgo  records  , indicated 

that  the  stress  duration  was  about  one  microsecond,  in  agreement  with 

tho  result  obtained  from  simple  calculations  using  the  measured  clastic 

wave  spoed  and  the  thicknoasea  of.  tho  specimen  and  flyer  plate.  Thus 

5 

crack  propagation  velocities  of  2  x  10  cm/soc  are  indicated — an 

interesting  result,  since  it  is  approximately  one-third  of  the  measured 

longitudinal  wave  velocity  for  novaculite,  and  is  thus  in  agreement  with 

theoretical  estimates  of  the  maximum  crack  velocities  of  brittle  i.  leri- 
8-11 

als.  Materials  such  as  Armeo  iron  and  beryllium,  investigated  in 

13,16 

other  projects,  fractured  in  a  brittle  manner  at  high  strain  rates, 

but  exhibited  viscous  crack  propagation  and  maximum  crack  speeds  well 

below  c./3. 
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VI  COMPARISON  OK  PREDICTED  AND  BCPF3U MENTAL  RESULTS 

A  tost  of  the  predlctlvo  capability  of  the  model  was  mode  by 
calculating  the  fragment  size  distribution  of  Experimont  53,  and  then 
comparing  tho  predicted  roault  with  the  experimentally  measured  result. 

Experimental  Conditions 

Experiments  52  and  53  were  carried  out  under  identical  conditions 
(Table  2).  Both  novsculite  specimens  wore  encased  in  aluminum  to  retain 
th©  fragments,  and  both  were  impacted  at  a  velocity  of  about  4ft  m/sec. 
Specimen  52  was  sectioned  afterwards  to  reveal  the  location  of  cracks 
and  fragments  (Figure  4d),  whereas  the-  aluminum  encasement  was  carefully 
removod  from  Specimen  53  and  the  size  distribution  of  the  fragmonts  was 
determined  (Figures  5  and  6). 

in  both  experiments  tho  novaculito  specimens  were  bended  on  the 
impact  side  to  tho  aluminum  with  epoxy;  tho  rear  surface  however  wan 
unbonded,  and  it  appears  therefore  that  a  gap  perhaps  10  to  30  urn  wide 
existed  there.  It  was  not  roaliz*  .<  at  the  time  of  the  experiment  how- 
critical  tho  gap  at  the  rear  interface  would  bo,  and  consequently  no 
stops  wore  taken  to  either  eliminate  or  measure  tho  gap.  Subsequent 
wave  propagation  calculations,  however,  showed  that  the  presence  of  a 
gap  has  a  largo  influence  on  tho  location  and  magnitude  of  tho  peak 
tensile  stresses  and,  thoroforo,  on  tho  fracture  behavior. 

Calculat lor.al  Conditions 

Computations  made  at  stress  levels  bolow  that  required  for  fracture 

for  the  situation  whore  a  bond  existed  at  the  rear  interface  predicted 

2 

a  peak  tonsion  of  about  160  MN/m  (23,200  psi),  which  first  appeared 

near  tho  rear  interface.  For  the  situation  of  an  unbonded  interface 

2 

with  no  gap,  a  peak  tension  of  100  MN/m  (14,500  psi)  occurred  first 
noar  the  Interface  and  dropped  off  rapidly  as  a  gap  formed.  For  tho 
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case  whore  an  initial  gap  existed  at.  the  interfaco,  the  roar  faco  of 

tho  8pocimon  behaved  as  a  free  surfaco,  and  tensile  stresses  of  about 
2 

100  MN/ra  (14,500  psl)  wore  predicted  throughout  tho  specimen  thickness, 
arising  first  at  , about  one-third  of  the  specimen  thickness  from  the 
rear  interface.  This  last  stress  distribution  is  in  accord  with  the  - 
cracking  patterns  observed  in  Figure  4d,  and  so  we  assumed  in  calculating 
the  fragment  size  distribution  tor  Experiment  53,  that  a  0,0025-cm  gap* 
ex* s ted  at  the  roar 'interface.  .  ...  .  '.•■'!•* v* ':"J' '•  ;l 

In  attempting  to  predict  the  fragment  size  distribution  for 
Experiment  53,  we  assumed  8-sidcd  fragments  and  fragment  sizes  equal  to 
tho  sizes  of  tho  cracks  at  the  time  of  coalescence,  i.e.,  wo  let 
8  b  1/4  and  y  =  1.0.  The  calculation  was  stopped  aftor  five  reverbera¬ 
tions  of  the  tensile  waves. 

Predicted  versus  Observed  Results 

The  results  of  the  calculations  are  presented  in  Figure  11,  which 

shows  that  fragmentation  occurred  at  three  positions  through  the  specimen 

thickness.  The  coalescence  criterion  was  choson  by  comparing  this 

computed  result  with  tho  observed  crack  patterns  of  Figure  4.  This 

comparison  indicates  that  tho  total  crack  range  volume  for  coalescence 
3 

should  be  about  3  cm  ,  whicn  implies  a  value  for  T  of  IT/3.  In  view 

c 

of  our  previous  estimates  of  the  T  parameter  in  Chapter  IV,  this  value 

■  c 

is  somewhat  low,  but  several  calculations  with  different  values  of  T’ 

c 

show  that  there  is  no  significant  effect  on  the  fragment  sizes,  -Rarefac¬ 
tion  waves  interacted  to  produce  tension  first  near  Cell  40,  and  tho 
stress  duration  was  sufficiently  long  that  enough  crack  growth  occurred 
for  the  total  crack  range  to  equal  the  cell  volume.  Hence  Cell  40 
fragmented.  During  the  third  tensile  reverberation,  Cell  44  fragmented 
and  during  the  fifth,  Cell  21  also  fragmented.  Cell  22,  which  exceeded 
the  coalescence  condition  by  nearly  as  much  as  Cell  21,  did  not  reach 
tho  f ragmontatlon  condition.  Thus  tho  calculations  predict  two  heavily 
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fractured  position#  in  tho  vicinity  of  each  facu  of  the  specimen — a 
prediction  that  is  In  reasonable  accord  with  tho  experimental  result  as 
shown  In  Figuro  4d. 

Although  the  fracture  pattern  In  Spocimen  52  (Figure  4d>  was  prodlctod 
well  by  the  code  calculation,  the  agreement  may  be  somewhat  fortuitous. 
Wheroas  tho  calculations  showed  cracking  and  fragmentation  occurring 
first  at  two  locations  near  the  roar  surface  and  later  at  locations 
nearer  the  impact  surface,  the  aeries  of  micrographs  of  Figure  4  Implies 
that  cracking  occurs  first  near  tho  impact  surface  and  later  near  the 
rear  surface.  Thus  if  the  micrographs  of  Figure  4  may  be  thought  of  as 
a  time  sequence,  the  experimental  observations  are  opposite  from  the 
fracture  behavior  predicted  by  tho  code. 

To  computo  tho  fragment  size  distribution,  the  fragments  in  all  the 
cells  wore  summed.  In  Cells  21,  40,  and  44,  tho  fragments  were  counted 
and  sized  as  described  in  Chapter  IV.  In  the  other  cells,  where  the 
total  crack  range  did  not  meot  the  fragmentation  critorlon  but  did  moot 
the  coalescence  criterion,  wo  computed  the  fragments  in  the  usual  we-  and 
disregarded  the  unfragmented  matorial .  Tho  computed  fragment  size  distri¬ 
bution  is  presented  in  Figure  12,  where  the  distribution  measured 

* 

experimentally  is.  also  given  for  comparison, 

Tho  calculated  and  experimental  curves  in  Figure  12  do  not  coincide; 
but  they  are  qualitatively  alike.  Agreement  concerning  the  sizes  of  the 
large  fragments  is  good,  and  the  curves  have  tho  same  shupo  for  Inter- 
medlate-slzod  fragments.  There  appear  to  bo  too  mony  computed  fragments 

*  We  attempted  to  obtain  an  extra  point  on  the  experimental  curve  by 
measuring  the  radii  of  tho  three  la;  »r.st  frogments.  This  point  is 
included  in  Figure  12  and  used  to  guide  the  extension  of  the  experi¬ 
mental  curve  (broken  line).  The  three  largest  fragments  weighed 
0.4367  g. 
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in  tho  rango  of  c  »  0.03  cm.  Thin  suggoatn  that  changes  should  bo  made 
in  either  tho  growth  or  fragmentation  process  for  the  smeller  size  creeks. 
Tho  cutoff  in  the  computed  distribution  curve  at  small  fragment  sizes 
exists  bocauso  no  fragments  wore  computed  for  those  sizes  since  inherent 
flaws  smaller  than  about  0.02  cm  wore  not  activated. 

No  attempt  was  made  to  repeat  tho  calculations  using  diffe  snt  values 
of  3  and  y.  However,  it  is  likely  that  bettor  agreement  could  bo  obtained, 
and  It  would  bo  valuable  in  future  work  to  study  the  consequences  of  vari¬ 
ations  In  those  parr  voters  for  cases  where  more  experimental  data  ore 
avallablo. 

The  above  fracture  and  fragmentation  model  applies  only  to  uniaxial 
strain  loading  conditions.  However,  tho  generalization  to  two  dimensional 
axially  symmetric  loading  goometri.es  such  as  those  obtaining  in  many 
blasting  and  drilling  situations  is  possible,  and  is  in  fact  currently 
in  progress  under-  Contract  DAAD05-73-C-0025  with  the  Ballistic  Research 
Laboratory,  Aberdeen,  Maryland, 

In  this  report  all  lnclpK'.t  flows  aro  assumed  to  lie  normal  to  the 
direction  of  wave  propagation.  However,  our  present  fracture  model  allows 
an  initial  flaw  orientation  distribution  which  is  maintained  during 
propagation.  In  this  case  the  driving  and  opening  stress  for  a  given 
crack  is  the  component  normal  to  the  crack.  This  option  was  not  exercised 
in  tho  present  work  because  tho  incipient  flaws  nearly  nil  lay  in  a  plane 
normal  to  the  direction  of  wove  propagation. 
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FIGURE  11  VARIATION  OF  COMPUTED  . \)TAL  CRACK  RANGE  WITH  POSITION  FOR  SPECIMEN  S3 
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APPENDIX  I  Cl W RACT ERI ZAT 1  ON  OK  OTHER  ROCKS 
■* 

Four  rock  typos  wore  used  in  thi  t  study.  Arkansas  novneulito  was 
tho  material  most  complotoly  investigated,  but  characterization  work 
was  also  performed  on  Sioux  (Jaipec)  quartzite,  Holoton  limestone  (Pink 
Tennessee  marble),  and  Westerly  granite.  This  section  describes  the 
microatructuros  and  defect  structures  of  the  four  rock  types  and  presents 
tho  results  of  measurements  of  their  physical  and  mechanical  properties. 


Microstructures 
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Figure  13  shows  a  polished  nurface  of  Arkansas  novacullte  that  was 
otchod  for  2  minutes  in  40%  HF  at  room  temperature  to  roveal  the  grain 
structure.  •  Black  areas  are  holes  whore  inherent  flaws  intorsoct  the 
surface  or  where  grains  have  been  removed  during  tho  polishing  process; 
bright. areas  are  causod  by  reflected  light  from  intornal  flaw  surfaces. 

By  exploiting  tho  translucency  of  novaculito  and  focusing  into  tho 
material  to  a  depth  of  about  100  ura,  we  found  that  the  inherent  flaws 
ox.'  t  predominately  in  two  shapes,  penny-shaped  and  pencil-shaped. 

Boat  results  were  obtained  by  viowing  the  specimens  in  reflected  polar¬ 
ized  light  through  a  microscope  slide  und  with  an  oil  film  of  matching 
refractive  index  (n  =  1.55)  on  the  specimen  surface.  Figure  14  shows  n 
penny-shaped  and  a  pencil-shaped  flaw  slightly  below  tho  surface  of 
polish  and  inclined  at  an  angle  to  it,  so  that  only  a  section  of  each  is 
in  focus.  The  planes  of  tho  rather  homogeneously  distributed  penny¬ 
shaped  flaws  are  roughly  parallel  to  one  anothor,  and  most  of  the  pencil¬ 
shaped  flaws  arc  inclined  at  about  45  degrees  to  theso  planes. 

One-inch  cubes  of  Sioux  quartzite,  Westerly  granite,  and  pink 
Tennessee  marble  were  cut  from  the  large  blocks  received  from  the 

*  Supplied  by  the  Property  Determination  Research  Support  Group,. 

Twin  Cities  Mining,  Research  Center,  Bureuu  of  Mino3,  Twin  Cities, 
Minnesota  55111. 
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'FIGURE  13  MICROST RUCl URE  OH  ARKANSAS  NOVACUl.l TG  SHOWING  THE  FQUIAXED 
OUART^  CHAINS 
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<»l  PHOTOGRAPH  WITH  POLARIZED  LIGHT  (b)  SCHEMATIC  DEPICTION 
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Bureau  of  Minos  and  polished  on  three  perpendicular  sides  in  preparation 
for  petrographic  -.•‘xamination,  Photomicrographs  allowing  the  grain 
structures  of  the  throe  rock  types  are  presented  in  Figure  15. 


■j' 

•? 


i 

f* 

r 

i 


The  Sioux  quartzite  is  relatively  pure,  dense,  and  homogeneous. 

Largo  cracks,  pores,  and  faults  are  noticeably  absent.  The  grains  are 
oquiaxed,  randomly  oriented,  and  about  30  times  larger  than  those  in 

Arkansas  novaculite  (average  grain  diameter  Is  of  the  order  of  300u). 

•  •  •'  \  , 
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As  indicated  by  the  pronounced  relief  of  polished  surfaces,  Westerly 
granite  consists  of  hard  grains  (quartz)  In  a  softer  matrix  (microcline 
and  plagioclase) .  The  quartz  groins  are  generally  irregular  with 
diameters  often  exceeding  1000u.  The  dark  blotite  phase  is  randomly 
oriented. 

Tho  grain  size  in  the  marble  ranged  from  very  small  (■"■'lOu)  to  very 
large  <3000u)  and  was  easily  discernible  in  3/4  polarized  light.  A 
large  majority  of  the  grains  exhibited  pronounced  twinning.  No  preferred 
grain  orientation  was  evident. 

Inherent  Flaw  Structures 

The  inherent  flaw  structure  of  Arkansas  novaculite  was  easily 
discernible,  and  it  was  possible  by  counting  and  measuring  flaw  traces 
on  polished  surfaces  to  dotermine  quantitatively  the  inherent  flaw  size 
distribution.  This  effort  was  described  fully  in  Chapter  V.  The 
inherent  flaw  structure  of  Sioux  quartzite,  Westerly  granite,  and  pink 
Tonnessoo.  marble,  however,  was  much  more  difficult  to  see.  Nearly  all 
flaws  in  these  rocks  are  associated  with  grain  boundaries  and  could  bo 
detoctod  only  by  focusing  painstakingly  up  and  down  with  the  optical 
microscope  at  magnifications  greater  than  100X.  Occasional  transgranular, 
cracks  were  observed  in  the  feldspar  grains  of  the  Westerly  granite. 
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FIGURE  15 


MICROSTRUCTURES  OF  (a)  SIOUX 
QUARTZITE,  (b)  WESTERLY  GRANITE 
ANO  (<;)  FINK  TENNESSEE  MARBLE 


Speciul  vi.ewi.ng  and  crack  decoration  techniques  were  tried  in 
attempting  to  obsorvo  tho  flaw  structure.  Phase  contrast  photography 
and  scanning  electron  microscopy  proved  ineffective;  likewise  swabbing 
polished  rock  surfaces  with  silvor  nitrate  and  vacuum  impregnation  with 
an  organic  fluorescing  agent  to  decorate  the  microcracks  was  of  llttlo 
use.  Thermal  grooving  was  not  attempted,  but  soems.  of  doubtful  value 
■inco  the  flaws  are  associated  almost  exclusively  with  grain  boundaries 
which  themselves  should  bo  attacked  by  tho  thermal  grooving  process,  . 

It  might  be  fruitful  to  attempt  to  relate  grain  size  or  some  other 
readily  observable  characterizing  parameter  of  tho  grain  boundaries  to 
the  number,  sizes,  and  shapos  of  crack-liko  defects  between  the  grains, 

17 

A  procedure  recently  reported  by  Drace  at  al.,  which  uses  ion 
thinning  to  reveal  cracks  in  Westerly  granite  and  Rutland  quartzite, 
appears  promising,  flowovor,  we  wore  not  able  to  try  it  in  this  work. 

Physical  and  Mechnnical  Properties 

The  density  of  each  rock  t ■ po  was  monsurod  by  an  immersion  technique 
and  the  longitudinal  and  transverse  sound  wave  velocities  wore  determined 
by  the  tlrao-in-flight  method.  The  results  are  given  in  Table  4.  The 
density  measurements  were  all  near  tho  theoretical  value  and  indicated 
that  porosities  were  loss  than  1%.  Tho  Young's  moduli  for  the  four 
rock  types  were  calculated  from  the  measured  densities  and  sound  speeds 
according  to  the  relationship 


|~(l+v)  (l-2v)*| 
L  "  (1-v)  J 


(24) 


where  a  value  of  0.25  was  assumed  for  Poisson's  ratio,  and  the  resulting 
values  are  included  in  the  table. 
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Quanl-Statlc  Strength  Measurements 


Tho  tenoilo  strength  of  novncullto  and  Sioux  quartzite  under 

quasi-static  loading  conditions  was  determined  using  tho  Shi  expanded 
18  19 

ring  tost  *  shown  in  Figuro  16.  In  this  tost  hydrostatic  pressure 
acts  radially  against  tho  inside  wall  of  a  cylindrical  speclmon  to 
create  a  uniform  tangontlol  tensile  stress  in  the  specimen  wall,  Non- 
axlal  stressos  caused  by  misalignment  and  localized  stress  concentrations, 
which  normally  arise  from  gripping  or  supporting  the  test  specimens, 
are  eliminated  in  this  method. 
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Ring-shaped 

specimens  were  diamond  ground  from  oversized  blanks 

to  tho  following  dimensions: 

Arkansas 

Sioux 

•  * 

Novaculite 

Quartzite 

I.D.  (mm) 

50.8 

81.4 

O.D.  (mm)  58. 4  87.8  to  91.2 

Height  (mm)  7.62  26.4 

-4  -1 

All  strength  measurements  wore  made  at  a  strain  rate  o.f  3  x  10  sec  . 


Measurements  on  novaculite  were  token  both  in  isirund  under  a 
-5 

vacuum  of  10  torr.  Some  materials  show  higher  strengths  under  vacuum, 

indicating  that  stross  corrosion  is  an  importent  factor.  No  significant 

difference  was  seen  for  novaculite.  Results  of  15  tests  in  air  Indicated 

2 

a  tensile  Strength  of  44.1  ±  3.2  MN/m  (6400  ±  440  psi)  and  *2  tests 

2 

under  vacuum  indicated  44.9  ±  3.0  MN/m  (6500  ±  480  psi).  Tho  extremes 

2 

measured  were  50,7  and  38.6  MN/m  ,  The  results  of  nine  tests  on  Sioux 

2 

quartzite  In  air  indicated  a  tensile  strength  of  13.3  ±  2.3  MN/m 

(2660  ±  330  psi).-  Average  values  of  the  quasi-static  tensile  strength 

of  Westerly  granite  and  pink  Tennessee  marble  have  been  reported  by 

2  2 

Wawcrsik  and  Brown  to  bo  10.8  MN/m  (1575  psi)  and  8.14  MN/m  (1185  psi), 
respectively. 
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It  la  noteworthy  that  the  quasi-static  tensile  strength  of  novneulite 

is  the  sane  as  tho  high  strain  rate  tensilo  strength  moasurod  In  the 

Impact  experiments  In  Chapter  II.  The  strain  rate  of  the  dynamic  tests 

was  ostirahted  from  tho  oscilloscope  records  of  the  yttorbium  stress 

gagos  in  the  instrumented  experiments.1'  The  rise  time  of  the  stress 

pulse  is  about  a  tenth  of  a  microsecond,  the  elastic  strain  ia  about 
-4 

.7  x  10  ,  hence  the  strain  rato  is • estimated  to  havo  been  approximately 

3-1 

7  x  10*  sec  ..  Since  this  is  more  than  seven  orders  of  magnitude  higher 
than  the  strain  rates  of  tho  quasi-static  ring  toots,  tho  negligible 
difference  in  strength  values  indicates  that  the  fracture  strength  ia 
strain-rate  insensitive,  at  least  in  this  range. 


Stress  Wave  Measurements 

Four  attempts  were  made  to  measure  tho  load. history  in  specimens 
of  Westerly  granite  that  were  impacted  in  tho  gas  gun.  The  experimental 
conditions  were  similar  to  those  described  in  Chapter  II  and  Reference  1. 

Small  grid- like  gage?  made  of  ytterbium  wore  cemented  to  the  back 

surfaces  of  the  rock  specimens  and  backed  by  n  plexiglass  plato.  Because 

of  tho  plnxoresistlve  nature 'of  ytterbium,  a  change  in  resistance  is 

produced  in  tho  material  when  traversed  by  a  compressive  wave.  This 

resistance  change  is  measured  and  recorded  by  on  oscilloscope  as  a 

change  in  voltage,  and  is  subsequently  converted  to  a  rocord  of  stress 

22 

versus  time  by  means  of  an  established  calibration.  The  tensile  stress 
history  In  tho  rock  specimen  is  then  calculated  from  the  measured 
compressive  stress  history  in  the  uluxigluss  and  the  relative  shock 
impedances  of  the  two  materials.  (See  the  Appendix  to  tho  First  Annual 
Report). 

Four  specimens  of  Westerly  grant: o,  10  cm  square  by  0.625  cm  thick, 
were  impacted  by  flat  plexiglass  flyer  plates  0.165  cm  thick  nt  velocities 
of  27.4,  123,  3G.9  and  41.0  m/sec.  The  Impact  velocities  were  chosen 
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such  that  stresses  in  excess  of  the  dynamic  tensile  strength  wore,  , 
obtained.  Tho  purpose  was  to  record  the  fracture  signal,  <*  second 
hump  on  tho  stress  gage  record  caused  by  impingement  on  the  gago  of 
recomprossion  waves  emanating  from  internal  surfaces  as  cracks  form 
and  grow.  In  oil  four  experiments,  tho  impact  too  nonplanar  to 
yield  onalyzablo  results. 

Figure  17  gives  the  measured  stress  history  at  the  rear  surface  of 

specimen  1.260  impacted  at  41  m/sec.  A  clear  fracture  signal  was 

recorded,  hut  nonplanarity  of  impact  resulted  in  a  triangular  rather 

than  a  flat-topped  initial  compressive  wave,  and  therefore  a  direct 

and  unambiguous  determination  of  the  dynamic  tensile  strength  of 

Westerly  granito  is  not  obtainable  from  this  record.  Hc-.aver,  the  peak 

compressive  stress  waa  calculated  for  this  experiment,  and  using  tho 

results  of  this  record  an  upper  limit  on  tho  dynamic  tensile  strength 

2 

of  Westerly  granite  was  determined  to  bo  47  MN/m  (6800  psi). 
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FIGURE  17  EXPERIMENTAL  STRESS  HISTORY  NEAR  THE  REAR  SURFACE  OF  A 
WESTERLY  GRANITE  SPECIMEN  MEASURED  WITH  AN  YTTCRBIUM 
STRESS  GAGE 
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